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ABSTRACT 

The diverse properties of broad-line quasars appear to follow a well-defined main sequence along which the 
optical Ee II strength increases. It has been suggested that this sequence is mainly driven by the Eddington 
ratio (L/LEdd) of the black hole (BH) accretion. Shen & Ho demonstrated with quasar clustering analysis that 
the average BH mass decreases with increasing Ee II strength when quasar luminosity is fixed, consistent with 
this suggestion. Here we perform an independent test by measuring the stellar velocity dispersion cr* (hence 
the BH mass via the M-cr* relation) from decomposed host spectra in low-redshift Sloan Digital Sky Survey 
quasars. We found that at fixed quasar luminosity, ct* systematically decreases with increasing Ee II strength, 
confirming that Eddington ratio increases with Ee II strength. We also found that at fixed luminosity and Ee II 
strength, there is little dependence of cr* on the broad H/3 EWHM. These new results reinforce the framework 
put forward by Shen & Ho that Eddington ratio and orientation govern most of the diversity seen in broad-line 
quasar properties. 
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1. INTRODUCTION 

The observed multi-wavelength properties of broad-line 
quasars display a remarkable regularity in the sense that these 
properties are covariant to a large extent, indicative of sim¬ 
ple underlying physical mechanisms that are ultimately con¬ 
nected to the BH accretion process. The dominant trend is 
known as Eigenvector 1 (EVl), a physical sequence along 
which many quasar properties correlate with the strength of 
the optical Ee II emissio n. Since the discovery of EVl by 
iBoroson & GreenI ([1993), it has been frequently suggested 
that EVl is mainly driven by the Eddington ratio (oc L/M), 
based on the EVl behaviors in terms of X-ray properties 
(|V^g^Uayjj_996ULaoi]^ 19971) . Civ properties (IWills et alJ 
1200(1 ISulentic et al.ll2000h . and later on more quantitative ar¬ 
guments based on the virial BH m ass estimates that utilize the 
width of the broa d emission lines (lLaodl20()0l lBorosonll2002l 
iDong et ^120111) . 

A practical concern in the above arguments is the poten¬ 
tial effect of orientation on the observed broad H/3 width, 
and hence the virial BH mass estimates. There is now am¬ 
ple evidence (both theoretical arguments and observations) 
that the broad-line region has a flattened geornetry ( e.g., 

IWills & Brownell986tlRunnoe et al.ll2013HGaskellll200^ and 

references therein). Hence it is important to disentangle phys¬ 
ical effects (such as due to accretion) from geometrical ef¬ 
fects (such as due t o orientation) whe n studying quasars in 
the context of EV1. IShen & H^ (12014) proposed that, based 
on several novel statistical tests with the largest low-redshift 
quasar sample to date, the diversity of quasar properties in the 
2-dimensional EV 1 plane defined by optical Ee II strength {x 
axis) and broad H/3 EWHM (y axis) is mainly determined by 
the Eddington ratio (traced by Ee II strength) and orientation 
(traced by broad H/3 EWHM) at fixed quasar luminosity. 
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The primary methodology used in IShen & Hoi (12014) to ar¬ 
gue that EV 1 is driven by Eddington ratio is quasar clustering, 
which can differentiat e low-mass and h igh-mass quasars in a 
statistical sense (e.g.. IShen et alJl2009h . without using virial 
BH mass estimates that are prone to or ientation and other 
systematic uncertainties (e.g.. lShenll2013h . An alternative ap¬ 
proach is to measure the stellar velocity dis persion (tr^,) of 
quasa r hosts, and use the M-a^ relation (e.g. lGebhardt et aT] 
lioOOt lEerrarese & Merritl 120001) to infer the BH mass in 
quasars independently from the clustering test. 

In this Letter we perform the latter exercise. Using a 
large spectroscopic qua sar sample drawn from the SDSS DR7 
(lAbazaiian et alJl2009l) . we were able to measure host ct, as 
functions of quasar properties (such as luminosity, H^ width, 
and Ee II strength) by performing host/quasar spectral decom¬ 
position. We found that the host velocity dispersions follow 
consist ent EV 1 patterns as expected from the empirical frame¬ 
work in[Shen^&[H3(l20l3)- In ^we describe the data and our 
spectral measurements; we present the main results in ^and 
summarize our findings in ® We adopt a flat ACDM cosmol¬ 
ogy with Do = 0.3 and Hq = 70kms“'Mpc“* throughout. Eor 
ease of discussion, we will use the terms quasar and Active 
Galactic Nucleus (AGN) interchangeably. 

2. DATA AND SPECTRAL MEASUREMENTS 

Our quasar sample consists of all objects spectroscopically 
classified as broad-lin e quasa rs (CLASS=' QSO') from SDSS 
DR7 (lAbazaiian et al.ll2009lF l with a redshift cut of 0.01 < 
z < 0.8 such that we have a sufficient wavelength coverage 
to measure the host tr*. Prior to the redshift cut, we have re¬ 
moved spectra from bad spectroscopic plates (PLATSN2=0), 
and we have checked the consistency between redshift val¬ 
ues given in SDSS “spSpec*” files and those measured by the 
idlspecid pipelin^B- Object spectra are discarded if both of 
the two redshifts are flagged as potentially problematic and 
inconsistent with each other. Our parent quasar sample con¬ 
tains 31346 unique objects. Pig. [T] displays the redshift and 

^ http://das.sdss.org/www/html/ 

® Publicly available at http://www.sdss3.org/svn/repo/idlspec2d/ 
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Figure 1. Top: Distribution of quasars in the luminosity-redshift plane. The 
continuum luminosity is estimated at restframe 5100 A from the full spectrum 
before host deco mposition). Objects included in the DR7 quasar catalog 
Shen et alJ201 111 are shown in red dots, while the black dots show additional 
quasars included in our sample. The better coverage of our sample in the 
low-redshift (low-luminosity) region is clearly shown here. Bottom: Redshift 
histograms of our sample (black) and the DR7 quasar catalog (red). 

luminosity distributi on of our sample. Compared to the SPSS 
DR7 quasar catalog (ISchneider et al.ll2010l:lshen et al.ll2Mlh . 
our sample expands the dynamic range in luminosity by in¬ 
cluding low-luminosity quasars that failed to meet the lumi¬ 
nosity threshold in the DR7 quasar catalog. The wavelength 
coverage of SDSS spectra is 3800-9200 A, with a spectral res¬ 
olution of R ~ 2000, and a pixel size of 10“^ in logjg A, which 
corresponds to 69kms“'. If there are duplicate spectra for 
the same quasar, the one with the higher S/N is adopted for 
spectral measurements. We do not distinguish radio-quiet and 
radio-loud quasars in this work, as the l atter population is only 
^ 10% of the e ntire quasar populat io n (IShen et al.l 201 Ih. 

We follow IVanden Berk et all (120061) and Shen et aP 
(l2008bh to decompose the SDSS spectra into pure-host 
and pure-quasar spectra, using eigenspectra of quasars and 
galaxies constructed using the P rincipal Compon ent Anal¬ 
ysis (PCA) upon SDSS spectra (lYip et akl l2004al[bl) . This 
method assumes that the composite (quasar-i-galaxy) spec¬ 
trum can be described by the combination of two indepen¬ 
dent sets of eigenspectra de rived from pure galaxy an d pure 
quasar samples. Empiricallv. lVanden Berk et al.1 (120061) found 
that a PCA decomposition using only the first few galaxy 
and quasar eigenspectra can reasonably recover the proper¬ 
ties of the host galaxy, provided there is a significant contri¬ 
bution from the host galaxy in the composite spectrum. The 
spectral decomposition is deemed successful if the host frac¬ 
tion, fn = /host/(/host + /qso), within 4160A < A < 4210A is 
greater than 0.1 and less than 0.9. The lower boundary is 
fou nd necessary in recoverin g the actual host galaxy proper¬ 
ties (IVanden Berk et al.l2006h . while the upper boundary is set 
to ensure successful measurements on quasar properties. We 
obtain the decomposed galaxy (quasar) spectrum by subtract¬ 
ing the reconstructed quasar (galaxy) spectrum from the orig¬ 


inal spectrum, and use t he original spectral error arrays in the 
spectra l measurements (IVanden Berk et akl 120061: 1 Shen et akl 
l20081j) . We obtained successful decompositions for 20254 
objects. 

We then fit the decomposed host galaxy spectrum with the 
publicly available code vdispfit included in idlspecld to 
estimate the host tr*. vdispfit is a template-based fit¬ 
ting method that performs the fits in the pixel space, where 
the templates consist of the first five PCA ei genspectra con- 
structed from the echelle stellar spectra in iMoultaka et al.l 
(120041) and cover rest-frame wavelengths from 4125 A to 
6796 A. During the fitting we mask narrow emission lines ex¬ 
cited both by the galaxy and by the quasar, as well as the spec¬ 
tral region around the quasar broad H/3 line (4760-5020 A), 
where the host-quasar decomposition often shows significant 
residuals due to the small number of quasar eigenspectra used. 
We allow a small wavelength shift (within ±500kms“') in 
the fit to account for the uncertainty in the systemic redshift 
estimate. A fifth-order polynomial is also added to account 
for the broad-band continuum shape. The fit is restricted 
to the rest-frame wavelength range of 4125-5350 A, which 
includes numerous stellar absorption features such as the G 
band (4304 A), the Mg Ib AA 5167,5173,5184 triplet, and Ee 
(5270 A), but excludes the Ca Hh-K AA 3969,3934 region. 
During the fit the templates are broadened to the native SDSS 
spectral resolution before the minimization, and the reported 
(7* values are already corrected for instrument broadening. It 
has been shown that cr, values measured by vdispfit are 
consistent with t hose measured by other fitting methods (e.g., 
IShen et aT]l2015l). We have also tried the empirical approach 
in IShen et al. (I200^bh for aperture correction in ct*, but found 
no significant changes in our results. Therefore we use the 
directly measured cr* in our following analysis. 

We restrict our analysis to objects with well measured cr* by 
vdispfit ((T*/A(T* > 3). Our final sample includes 9999 
quasars with meaningful cr* measurements. We have com- 

J iared o ur cr* measurements with those reported in IShen et alJ 
20081j) for the common objects, and found consistent results. 
The success of decomposing the host spectrum and measur¬ 
ing (7* depends strongly on the host-quasar contrast, which 
in turn introduces a luminosity-dependent selection bias such 
that more luminous quasars are less likely to yield a successful 
(7* measurement. We will return to this point in ^ The typ¬ 
ical S/N per SDSS pixel for objects with successful cr* mea¬ 
surements is ^ 6 in the decomposed host spectrum, and as a 
result many of our cr* measurements have large uncertainties 
(as much as ~ 30%, possibly even larger given the system¬ 
atic uncertainties). Nevertheless, the average cr* for binned 
objects should be well determined. Throughout this work we 
only use the average (estimated by the median) values of cr, 
in our binned analysis instead of individual measurements. 

In addition to cr*, we also measure the quasar spectral prop¬ 
erties in the region, including the continuum luminosity 
T5100.AGN = ALa at restframe 5100 A, [O III], Ee II and broad 
p roperties, using the metho dology outlined in earlier work 
(e.g.. IShen et alJl2008ai 120111) . These quantities are neces¬ 
sary to investigate the trends of cr* in terms of the EV1 cor¬ 
relations. Among those objects with meaningful cr* results, 
7922 have good quasar measurements (e.g., the broad H/3 line 
is detected at > 3a), which form the basis of our following 
analysis. We have visually inspected the objects that failed 
to produce good broad-line measurements but have good cr* 
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Figure 2. Median values of host-galaxy ct* in the EV 1 plane. Columns cor¬ 
respond to increasing Speii from left to right, and rows correspond to increas¬ 
ing broad H/3 FWHM from bottom to top. The total number of objects in 
each Rpeii—FWHM h /3 bin is marked in each panel. The missing data points 
con'espond to bins containing less than 50 objects. At fixed luminosity, there 
is a general decreasing trend of it* with increasing Speii, suggesting that the 
Eddington ratio (L/M) increases with Speii. On the other hand, there is little 
dependence of cr* on the broad H/f FWHM at fixed /?pcii and luminosity. 

measurements, and found that ^ 60% are due to pipeline mis- 
classifications (the object is a galaxy or a type 2 AGN) and 
the remaining are due to insufficient S/N or weak broad H/3 in 
the fitting region (many of them show broad Ha). The mea¬ 
surements for our final sample are provided in a supplemental 
online EITS table. 

As a sanity check, our measurements reproduce the 
well known EVl relations reported in previous work 
(iBoroson & Greenlll992l : IShen & Holl2014ft . and extend these 
results to fainter quasar luminosities enabled by our sample. 

3. RESULTS 

We show our main result in Eig. |2] Our quasar sample with 
cr* and broad-line measurements was divided into subsamples 
in the two-dimensional EV 1 plane of Ee II strength, measured 
by the ratio of Eell EW within 4434^684 A to broad H/? 
EW (Rpeii = EWpeii/EWH/s), and the EWHM of the broad 
H/3 (EWHMh/j). In each Rpeu-EWHMn^ bin, we further 
divide the sample into different quasar luminosity bins. Sta¬ 
tistically speaking, higher luminosity quasars are on average 
more massive BHs, hence we need to compare the cr* values 
at fixed quasar luminosity. 

Eig.|2]demonstrates the behavior of cr* along the two axes of 
the 2D EV 1 plane. In each EWHMh /3 row and at fixed quasar 
luminosity, there is a systematic trend in the average cr* as 
Rpeii increases; quasars with high Rpeii generally have lower 
host stellar velocity dispersiorQ This suggests that the BH 

^ IGreene & Hd <2006f) have shown that the cr* measurements could be 
biased in quasar's with small-width broad Fe II emission when the fitting is 
restricted to certain spectral windows (typically several hundred A), partic¬ 
ularly for high Eddington ratio objects (i.e., the high Rpcll end of the EVl 
sequence). Our cr* fitting is over a broad range of wavelengths, and the 
Fe II emission is largely removed with our host-quasar decomposition, both 


Figure 3. The median (with eiTor) value of cr* as a function of RfcII, for each 
luminosity bin shown in Fig. [2] Note that objects with different FWHMh /3 
are merged in those L-Rpeii bins. A general trend of decreasing cr* with 
increasing Rpdi is seen at all luminosities. The relatively less prominence in 
the trend at low luminosities is likely due to the selection bias against small 
O'* in our sample. 

mass systematically decreases with Rpeii when luminosity is 
fixed, hence the Eddin^on ratio increases with Rpep. Another 
notable feature in Eig.^lis that there is little dependence of cr* 
on EWHMh /3 when luminosity and Rpeii are fixed, while the 
standard virial BH mass estimates predicts a strong vertical 
segregation in BH mass. The lat ter observation is c onsistent 
with the framework suggested bv IShen & Hoi (120141) : the dis¬ 
persion in EWHMh^ at fixed luminosity and Rpeii is largely 
due to orientation effects in a flattened BLR geometry, rather 
than intrinsic dispersion in BH mass. Therefore cr* remains 
more or less constant in different EWHMH,a bins when lumi¬ 
nosity and Rpeii are fixed. Given the insensitivity of cr* on 
EWHMh /3 in the 2D EV 1 plane, in the following analysis we 
only consider trends of cr* with luminosity and Rpeii by merg¬ 
ing objects with different EWHMh /3 in each L-Rpeii bin. 

Eig. [3] summarizes the dependence of host cr* as a function 
of Rpeii at different quasar luminosities. Again, the general 
trend of decreasing average cr* with increasing Rpeii is seen 
at all luminosities. However, one caveat is that the criterion 
to select well measured cr* (see ® and the SDSS spectral 
resolution may induce a bias against low cr* values. This may 
be responsible for the less prominent cr* -Rpeii trend seen for 
the low-luminosity quasars (with low cr* values) in Eig. [3 

We point out that at fixed quasar luminosity, the fraction of 
quasars with measurable cr* does not show significant varia¬ 
tion in each Rpeii-EWHMn/s bin, hence there is no selection 
bias in the observed trend of cr* with Rpeii or EWHMH,a- On 
the other hand, cr* is generally more difficult to measure in 
higher-luminosity quasars, and those objects with measured 
cr* may be biased towards massive hosts (hence higher cr*). 

of which help remed y for this caveat. Even if we were fitting cr* similarly as 
IGreene & Hd 420061) . for th e ranges of Eddington ratios and cr* relevant in 
this work, the simulations in IGreene & Hd 420061) suggest that any potential 
bias in cr* measurements is on the < 10% level, which is far too small to 
account for the trends in cr* seen in Figs.[2]an d[3] 
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In other words, the apparent trend of increasing cr* with lumi¬ 
nosity may be subject to this luminosity-dependent selection 
bias in terms of cr* measurement. 

To assess the importance of this potential luminosity 
bias in ct* measurements, we make simple estimations us- 
ing the observed M - relation in local inact ive galaxies 
(iKormendv & Holl2013t [McConnell & Mall2M3h . Assuming 
a constant Eddington ratio A = Tboi/TEdd in each /?Feii bin 
(where Lead = 1.26 x 1O^*M/M0ergs“' is the Eddington lu¬ 
minosity), we have an expected relation: 


log( 


CT ^ 


200km s 


, , E5100.AGN 1 0-7 1 n \\ 

log --37.1-a- logA) 

' ergs ^ ' 


( 1 ) 

where we have assumed logLboi = IOLsioo.agn, and converted 
BH mass to cr* using the local M- cr* relation: 


with a = 8.49 and /3 = 4.377 a s in IKor mendv & Hoi (120131) . 
and a = 8.32 and (3 = 5.64 as in iMcConnell & Mai(l20 1 3h ' 
The predications are shown in Eig. |4] and compared with 
data. We simply adjusted the value of the average Edding¬ 
ton ratio (log A) in each /?Feii panel to best match the predic¬ 
tions with the data (as judged by eye). The data points for the 
highest luminosity bins significantly deviate from the predi¬ 
cations, suggesting there might indeed be a strong selection 
effect in the highest luminosity bins that leads to overesti¬ 
mated cr* on average. In addition, as we mentioned earlier, 
the tr* measurements at the low-luminosity (1ow-ct*) end may 
also be biased high due to our measurement quality cut and 
the SDSS spectral resolution, which will flatten the relations 
shown in Eig. |4] Nevertheless, the general agreement in the 
trends between data and simple model predictions suggests 
that higher-luminosity quasars likely have intrinsically larger 
cr* (hence more massive BHs) on average when /?Feii is fixed. 
This is consistent with the idea that /?Feii is a good proxy for 
the Eddington ratio. As expected, the normalization of (log A) 
required to match the predictions with data increases as 7?Feii 
increases; these Eddington ratio estimates, while subject to a 
number of caveats, are also reasonable values for the low-z 
QlltiisiAi S ill our sample. In a companion paper (Sun & Shen, in 
preparation), we will explore the luminosity dependence of 
more carefully with coadded spectra in the EV1 plane, which 
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Figure 4. The median value of a* as a function of luminosity in different 
ranges of f?Feii (different panels). Objects with different FWHMh^ in a given 
L-)?Feii bin are combined since there is little dependence of tj* on FWHMh/? 
when luminosity and ^Feii are fixed (Fig.[2)- In addition to our data points, 
we show simple relations expected from the observed local M-a* relation, 
assuming a constant Eddington ratio (logA) in each i?Feii bin (e.g., Eqn.^. 
The dotted and das hed lines represen t the r esult s using the observed M — 
cr+ relation given by IKormendv & j2Q13l) and IMcConnell & Mai <20131) . 
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will allow us to better measure the average host velocity dis¬ 
persion as a function of quasar luminosity that is less suscep¬ 
tible to the incompleteness of individual measurements. 

4. DISCUSSION AND CONCLUSIONS 

We have studied the host stellar velocity dispersion in 
quasars as functions of quasar properties in the context of 
EV1, using a spectral PCA decomposition technique to sep¬ 
arate the host galaxy and quasar components in integrated 
SDSS spectra. Our sample includes all low-redshift (z < 0.8) 
quasars selected from the SDSS DR7, and extends to fainter 
luminosities than thos e included in the DR7 quasar catalog 
([Schneider et al.ll2010h . 

The trends of host stellar velocity dispersion tr* with EV 1 
properties revealed in this study provided new insights on 
the origin of the diversity of quasars in EVl space. We 
found that at fixed quasar luminosity, cr* on average de¬ 
creases with the optical Ee II strength, suggesting the average 
BH mass (Eddington ratio) decreases (increases) with Ee II 
strength, confirming the early sugg estion that Eddington ra¬ 
tio drives the EVl sequence (e.g., iBoroson & Green! [19921: 
lBorosonll2002l) . On the other hand, at fixed quasar luminos¬ 
ity and Ee II strength, there is little dependence of on the 
broad H/g EWHM , which is consistent with the suggestion in 
IShen & (120141) that orientation plays a significant role in 

the dispersion in H/3 EWHM at fixed Ee II strength. Although 
a quantitative mapping between Ee II strength and Eddington 
ratio is beyond the scope of this paper (which requires a more 
careful treatment of the selection effects and scrutiny on sev¬ 
eral additional assumptions), these finding s on provide fur- 
ther support to the empirical framework in IShen & Hoi (120141) 
that Eddington ratio and orientation govern most of the diver¬ 
sity in quasar properties at fixed luminosity. 

Establishing Eddington ratio as the physical driver for the 
many covariant properties of quasars has important implica¬ 
tions on the accretion physics of quasars. The change in ac¬ 
cretion rate may modify the st ructure and the o utput spec¬ 
trum of the accretion flow (e.g.. lWang et al.ll2014^ . and mod¬ 
ulate the photoionization processes both in the nuclear region 
and in the more distant narrow line region, and therefore lead 
to the coherent patterns we observe as the main sequence of 
quasars. The empirical evidence presented here and in earlier 
work well motivates a theoretical investigation on the connec¬ 
tion between quasar phenomenology and accretion processes. 
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SDSS Web Site is http://www.sdss.org/. 
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